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Abstract 

Here we continue to elaborate properties of the relativistic mean-field based model 
(SHMC) proposed in ref. [6] where hadron masses and coupling constants depend 
on the cr-meson field. The validity of approximations used in [6] is discussed. We 
additionally incorporate contribution of meson excitations to the equations of mo- 
tion. We also estimate the effects of the particle width. It is demonstrated that the 
inclusion of the baryon-baryon hole and baryon-antibaryon loop terms, if performed 
perturbatively, destroys the consistency of the model. 



1 Introduction 



In recent years there has been a great interest in the description of hadronic 
properties of strongly interacting matter. It is based on the fact that var- 
ious experiments indicate modifications of hadron masses and widths in the 
medium (see for example [1]). As expected previously, these changes are possi- 
bly related to a partial chiral symmetry restoration in hot and / or dense nuclear 
matter, cf. [2J. Later on it was realized that the connection between the chiral 
condensate of QCD and hadronic spectral functions was not as direct as it was 
originally envisaged. Nevertheless, the study of an in-medium modification of 
hadrons is an essential point of scientific programs at new heavy ion facili- 
ties at FAIR (Darmstadt) [3], NICA (Dubna) [I] and low-energy campaign at 
RHIC (Brookhaven) [5]. 
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Theoretical predictions for critical baryon density and temperature of the 
hadron-quark phase transition depend sensitively on the Equation of State 
(EoS) at high densities and temperatures. In [5] and here we focus on the 
study of the EoS of the hadronic matter. Any EoS of hadronic matter should 
satisfy experimental information extracted from the description of global char- 
acteristics of atomic nuclei such as the saturation density, the binding energy 
per particle, the compressibility, the asymmetry energy and some other. Def- 
inite constraints on hadronic models of EoS are coming from the analysis of 
direct and elliptic flows in Heavy-Ion Collisions (HIC). In addition to these 
constraints astrophysical bounds on the high-density behavior of /3-equilibrium 
neutron star matter should be applied, see [7J. 

In [6], we constructed a phenomenological Relativistic Mean-Field (RMF) 
based model that allows one to calculate particle in-medium properties and 
the EoS of hadronic matter in a broad density-temperature region. The valid- 
ity of this model was demonstrated for the description of heavy ion collisions 
in a broad collision energy range. Microscopically based approaches, as the 
Dirac-Brueckner-Hartree-Fock method, see [8], are very promising but need 
rather involved calculations. The model of ref. [B] is a generalization to finite 
temperatures of the RMF model developed in [§] and applied in [7J (KVOR 
model) for describing neutron star properties. 

Following ref. [9] we assume relevance of the (partial) chiral symmetry restora- 
tion at high baryon densities and/or temperatures [TD] manifesting in the 
form of the Brown-Rho scaling hypothesis [TT] : Masses and coupling constants 
of all hadrons decrease with the density increase approximately in the same 
way. In [6J we followed the simplest form of the scaling hypothesis and scaled 
the quadratic (mass) terms of a, u, and p fields, as well as the nucleon mass, 
by a universal scaling function <3> which was assumed to be dependent on the 
a mean field. In order to obtain a reasonable EoS, the meson- nucleon cou- 
pling constants were also scaled with the a mean field treated as an order 
parameter. Differences in scaling functions for the effective masses of the u- 
and p-fields and their couplings to a nucleon allowed us to get an appropriate 
density-dependent behavior of both the total energy and the nuclear asym- 
metry energy, in agreement with the constrains obtained from neutron star 
measurements, cf. 012] • 

Note that the idea of the dropping of the meson effective masses continues to 
be " a hot point" being extensively discussing in the literature. There exist 
works which simulate different modifications of the simplest form of the scaling 
trying to find an optimal ansatz. E.g., the model [32] introduces a common 
dropping of the uj, p effective masses, whereas a is treated differently, as purely 
classical field, i.e. the static space-independent order parameter. Scalings of 
the u, p effective masses on the one hand and the nucleon effective mass on 
the other hand are assumed to be different. Couplings are evaluated following 



2 



quark counting. As in [6] and in the given paper, u, p mesons are assumed 
to be coupled only to the classical a field, since in the quark model they are 
made of a quark and an antiquark, which couple oppositely to the vector 
field. A support for the common dropping of the N,a,u, p masses comes from 
lattice QCD in the strong coupling limit [H] where it was found that meson 
masses are approximately proportional to the equilibrium value of the chiral 
condensate. 

There exist models, which do not accept the idea of the dropping of the ef- 
fective meson masses at all. E.g., most of the RMF models continue to use 
the constant a, u, p effective masses. Some models introduce field interaction 
terms leading to an increase of the a, u, p effective masses with the increase of 
the nucleon density, e.g., see [T5][T6] . Ref. [TF] suggests an increase rather than 
a decrease of the p meson mass with increase of the temperature, motivating 
it by mixing of vector and axial mesons at finite temperature, that authors 
consider as an indication towards chiral symmetry restoration. Another mod- 
els simulate only the p width rather than a modification of the mass, although 
from general point of view a modification of the imaginary part of the self- 
energy of the resonance should stimulate a modification of the real part of the 
self-energy (effective mass), as a consequence of the Kramers-Kronig relation. 

In present paper, as in our previous paper [6] , we avoid discussion of these in- 
teresting theoretical questions. Instead we will follow the Brown-Rho scaling 
hypothesis in its simplest form confronting further the results of the model 
with the HIC data. Besides the nucleon and meson a, u and p mean fields, 
we included low-lying non-strange and strange baryon resonances, meson ex- 
citations <r(600), p(770), u(782) constructed on the ground of mean fields, 
and the (quasi) Goldstone excitations 7r(138), K(495), 7/(547) as well as their 
high mass partners in the SU(3) multiplet if* (892), 7/(958), and y?(1020). All 
corresponding antiparticles are also comprised. Interactions with mean fields 
are incorporated as well. In ref. [TS] it was shown that it is possible to re- 
produce particle scattering data when the lowest baryon octet and decouplet 
are assumed to be the only relevant degrees of freedom. Therefore we do not 
consider higher resonances within our model. 

In order to construct a practical model in [5|, we used several simplifica- 
tions. First, we assumed the validity of the quasiparticle approximation for 
all baryons and mesons. Second, we supposed that baryons and meson exci- 
tations interact only via a, uj and p mean fields. Thus, the fermion-fermion 
hole and the fermion-antifermion loop diagrams for boson propagators and 
the boson-fermion loop diagrams for fermion propagators were disregarded. 
Also, meson-meson excitation interactions were neglected. Thus, effectively 
excitations were considered as an ideal gas of quasiparticles. Treating meson 
excitations perturbatively we have omitted their contribution in the equations 
of motion. With this Scaled Hadron Mass- Coupling (SHMC) model we con- 
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structed the EoS as a function of the temperature and the baryon density and 
used this EoS in a broad density-temperature region to describe properties of 
hot and dense matter in heavy ion collisions. 

Note that the standard RMF models generalized to finite temperatures have 
been studied in the literature, e.g., see [T9f20f21j . In [19] temperature depen- 
dence was included only into nucleon distributions. A general treatment of 
meson excitations has been considered within the imaginary time formalism 
[2"U] and a more convenient real time formulation [2T]. We incorporate fluc- 
tuative terms expanding fields near their mean-field values. Simplifying we 
retain only quadratic fluctuations. Thus in our model the gas of excitations 
interacts only through mean fields. Within this approximation our results can 
be reproduced using above mentioned finite temperature quantum field theory 
techniques. 

In the present paper, we check the validity of different approximations as- 
sumed in [6] and consider several possibilities how the model can further be 
improved. In sect. [2] we introduce the SHMC model of [6]. In Sect. [3], the pres- 
sure functional of the model is constructed and the equations of motion are 
derived. Boson excitation terms are incorporated in the equations of motion 
and a comparison is made with the perturbative treatment carried out in [6]. 
Section [4] estimates the effects of finite particle widths. In Appendix A, we 
discuss differences in two possible treatments of the a meson field, first, as 
an order parameter (as in [B]) and second, as an independent variable, i.e., 
considering a on equal footing with other field variables (u> and p). Appendix 
B demonstrates problems which arise if the baryon loop terms are included. 
Fermion loop effects on the boson excitation masses are evaluated within a 
perturbation theory approach and arguments are given why these effects are 
not included into the SHMC model. 

In reality the nucleon self-energies have a momentum dependence which is not 
so small. It manifests itself in high energy heavy-ion collisions [22] and affects 
different properties of atomic nuclei [23]. P-wave pion- and kaon-baryon inter- 
actions may significantly affect properties of the pion and kaon sub-systems, 
see [2"4l25ll26j and refs. therein. As in [5] and in most of RMF models, here we 
continue to disregard the p-wave effects. 

A number of other important effects is not incorporated into our model. How- 
ever the full theoretical quantum field description of many strongly interacting 
hadron species can't be constructed in any case. Using RMF based models and 
their generalizations one should always balance between a realistic and practi- 
cally tractable descriptions. Thus we postpone with further generalizations of 
the SHMC model. Further improvements of the model will be done after it will 
pass the check in actual hydrodynamical calculations of heavy ion collisions 
in a broad energy regime, that is our future program. 
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2 About the SHMC model 



Following [H] we use the cr-field dependent effective masses of baryons 

m* b /m b = § b {Xa<?) = 1 - 9ab Xa cr/m b , be {b} (1) 
with the baryon set 

{b} = JV(938), A(1232), A(1116), E(1193), S(1318), E*(1385), S*(1530), and 
f2(1672) + all antibaryons. 

The mass terms of the mean fields are 

m* m /m m = |$ m (x CT cr)| , m E {m} = a,u,p, (2) 
where g ab are the crfr-coupling constants. 

For the sake of simplicity we scale all couplings g ab by a single scaling func- 
tion Xo-( cr )> an d all g^b, g P b by Xu{°) an d Xp{ a ) scaling functions, respectively. 
Therefore, all scaling functions depend only on a 0. The idea behind that 
is as follows. The cr-field can be interpreted as an effective field simulating a 
response of the wei-quark condensate. The change of effective hadron masses 
and couplings is associated, namely, with a modification of the quark conden- 
sate in matter. Thus, we consider the cr-field as a composite field, like an order 
parameter, whereas other meson fields are treated as fundamental fields. The 
cr excitations are then interpreted as fluctuations around the mean value of the 
order parameter. Similarly, long-scale fluctuations are treated in the Landau 
phenomenological theory of phase transitions. 

To single out quasiparticles (excitations) from the mean fields, one should do 
the following replacements in the Lagrangian: ujq = Uq + u/, Ro = Rq + Rq, 
uj = uj and p = p . Here u^, -Rq 1 are the mean (classical) field variables and 
w u> (Po) M ' (p'±Y are responsible for new excitations, Rq = p$. In [6] we con- 
structed a thermodynamic potential that besides mean-field terms includes 
the contribution of u and p excitations. By varying with respect to the fields 
a;^ 1 , Rq we obtain equations of motion from where the u)q(ct), Rq{&) fields are 
extracted and put back into the thermodynamic potential. A similar proce- 
dure has been used in a number of works, e.g. in [13]. Then, in contrast with 
[T3] . supposing a = cr cl + a', we expand the thus obtained effective potential 
in a' up to squared terms (contribution of o' fluctuations) and, varying the 
thermodynamic potential in cr cl , derive the equation of motion for the resulting 
order parameter. 

On the other hand, if the a field was treated on equal footing with u and 
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Ro, as it was done in the standard Walecka model, we would consider all 
three fields as independent variables. The comparison between two choices is 
performed in Appendix A. 

The dimensionless scaling functions and $ m , as well as the coupling scaling 
functions Xm, depend on the scalar field in the combination Xa( a ) a - Therefore, 
we introduce the variable 



/ = 9aN Xa °/m N . (3) 

Following [9] we assume approximate validity of the Brown- Rho scaling ansatz 
in the simplest form 



$ = = = = $ P = 1 - /, (4) 

using x<r — Thereby, in terms of a one obtains $(cr) = [1 + ff ff wff/mw] _1 . 
One could partially break the scaling, if it were required from comparison with 
the data. 



We keep the standard expression for the nonlinear self-interaction (potential 
U) of the RMF models, but now it is expressed in terms of the new variable 
/. Using ([3]) the potential U can be rewritten as follows: 



rr_ 4 / h f 3 i c f 4\ _ bm N (g aN Xo of c(g aN Xav) A /_v 
u —m N {- j + - / ) — - i - . {0) 

The presence of two additional parameters, "6" and "c", allows one to accom- 
modate realistic values of the nuclear compressibility and the effective nucleon 
mass at the saturation density. Extra attention should be paid to the fact that 
the coefficient "c" must be positive to deal with the stable ground state. Values 
of the parameters used in our SHM C model can be found in [6] . In Fig. [1] we 
present the dependence of nucleon (cf. Fig. 1 of [6]) and antinucleon optical po- 
tentials on the single-particle energy. Comparison is presented with predictions 
of the standard Walecka model (with only a and to mean fields). As it is seen, 
our model describes the nucleon optical potential in an optional way, better 
than the standard Walecka model. Differences in predictions of those models 
for antinucleon optical potentials are drastic. A phenomenological value of an 
antiproton optical potential is limited within the range —100 -. — 350 MeV 
|28j . in favor of the given model compared to the standard Walecka model. 
Predictions for antiprotons are very important in a light of future experiments 
at FAIR. 

There are mean-field solutions of the baryon and (T,u,p meson Lagrangian 
J2be{b} £-b + Em.e{m)^ F [§]• To these terms we add the Lagrangian density 
for all meson excitations 
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Fig. 1. Energy dependence of the nucleon (left) and antinucleon (right) optical 
potentials. Solid lines - predictions of our model and dash lines, of the original 
Walecka model. Shaded area shows uncertainties in extrapolation from finite nuclei 
to cold nuclear matter |27j . 

£ cx = E 4*, {ex} = 7r ± ' (138);ir ± ' ,i?°(495);r / (547); (6) 

exd{ex} 

ir ±,0 (892), t/(958), 0(1020); a\ uj\ p' . 

The set {g} = (tt, K, rj) is often treated as (quasi) Goldstone (index " g") bosons 
within the chiral SU(3) symmetrical models. Therefore, one may not scale their 
masses and couplings, as we have carried out for the mean fields a, ou, p, cf. 
the set A for couplings in Fig. 8 (left) of ref. [6] . On the other hand, one may 
observe, cf. [6.9] , that for the case of spatially homogeneous system the equa- 
tions for mean fields and thus their mean-field solutions do not change if one 
replaces the a, ujq, Rq fields by the scaled fields Xa°~> Xu u o an d XpRo, provided 
&b = &m = Xm, and x' p = X 2 P (x' P is the scaling function of the p — p interaction 
g p , see [6]). If one wishes to extend this symmetry to the case when Goldstones 
are included, in addition to the scaling of masses one should scale couplings, 
9mg = 9mg Xm, cf. set B in Fig. 8 (right) of [6]. In [6], we tested both possibil- 
ities g* mg = g mg and g^ g = g mg Xm, and referred to them as versions without 
and with scaling, respectively. We include interaction of (quasi) Goldstones 
with mean fields (for K and 77, for n it is small). As the result of this inter- 
action, at sufficiently large (overcritical) baryon densities there may appear 
mean field solutions for (quasi) Goldstone fields signaling of condensations of 
these fields. Values of critical densities are higher for the set B. Since there are 
no experimental indications of condensation of (quasi) Goldstone bosons in the 
heavy ion collision regimes, comparing our results with experimental data, as 
in [6], we will focus on the set B, where condensates do not occur. K* ±,Q , 77', 
are assumed not to couple with mean fields, since there is no experimental 
information for such a coupling. 

When the total Lagrangian is constructed, one can derive the equations of 
motion for every field. Even for low baryon density, the equations of motion 
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for a, uj and p allow mean-field solutions ctq, ujq, pjj. Therefore, we use 



cr = cr ; = <Lo; = -Ro ^a3 <W ( 7 ) 

We assume that the system volume is sufficiently large and surface effects 
may be disregarded. Thus, only spatially homogeneous RMF solutions of the 
equations of motion are considered. 



3 Improved description of meson excitations 



The thermodynamic potential density f2, pressure P, free energy density F, 
energy density E and entropy density S are related as 

E = F + TS, F[f,uj 0j R }=J2^ + n, n = -P, (8) 

i 

dF 

Summation index % runs over all particle species; are particle densities. 
Chemical potentials /i, enter into the Green functions in the standard gauge 
combinations gj + /ij. 

Thermodynamic quantities (JSJ) can be found from the energy-momentum ten- 
sor T^y which is defined by our Lagrangian. The energy density E and pressure 
P are given by the diagonal terms of this tensor 



E = (T 00 ) , P = \ (T u ) ■ (10) 

In [6], the energy was chosen as a generating functional. Here we will use the 
pressure functional since it is more suitable to treat meson excitation effects 
in the presence of the mean fields and baryon-loop contributions. 



3. 1 Pressure at finite density and temperature 



The pressure can be presented as the sum of the mean a- } uj-, p-field terms as 
well as of contributions of baryons and of all meson excitations. So we have 

P[f,u ,Ro}= £ P™ F [f,u ,Ro}+ W,u; ,R } 

+ -Pbos.ex. [/, • (11) 
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The first two sums are included in every RMF model but with a smaller set 
{b}, whereas the boson excitation term Pbos.cx. is constructed in [6] beyond 
the scope of the RMF approximation and will be further elaborated here. 

Although in our treatment of the a variable all terms in fill I) are functions only 
of / and T, we also present them as functions of Uo and Rq in such a way that 
values of the uJo(f) and Ro(f) mean fields can be found by minimization of the 
pressure. Then o>o(/) and Ro(f) are plugged back in the pressure functional 
that becomes a function of / only. The equilibrium value of / can be found 
by subsequent minimization of the resulting pressure in this field. 

In a self-consistent treatment, equations of motion for the mean fields render 



d 
duj 



P[f,<M=Q 



d 
dR 



P[f,Ro 



;i2) 



and 



±P[f^ (f),Mf)] = -HjP\fMf),Mf)] = ( 13 ) 

with pressure P given by eq. (JTTJ) . Since Pbos.ex. [/, ^o, Po] depends on the mean 
fields, its minimization produces extra terms in the equations of motion for 
the mean fields. In differentiating in (fT3|) we used ({12]) . This self-consistency of 
the scheme allows us to be sure of thermodynamic consistency of the model. 

In [6], excitations were treated perturbatively. Accordingly, we assumed that 
Pbos.ex. = Pbos.cx. [/ MP ,u, MF ,i2 MF ], where / M V MF ,P MF are found by mini- 
mization of the pressure without inclusion of the boson excitation term. Thus, 
equations of motion for mean fields that we used in [6] are: 



d 

d 
dRo 



£ P% F \f,u ,R ]+ E P b [f,u ,Ro] 
P^ F [f^ ,R }+ J2 Pb[f,u ,Ro\ 

m£{m} be{b} 



(14) 



and 



d_ 

df 



P^ F [f,oo ,Ro]+ E Pb[f,"o,Ro} 

mS{m} b£{b} 



d_ 

df 



E Pm F [f, UJ 0,Ro}+ E Pb[f,U ,Ro\ 
m£{m} be{b} 



0. 



(15) 
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Equation (fl4l) was used in differentiating in (|T5|) . Below in Figs.[2H5]we demon- 
strate how effects of a nonperturbative treatment of boson excitations, incor- 
porated in (1121) and f[T3"j) (a self-consistent analysis) and neglected in (114)) . 
(USD, affect results of the SHMC model. 

Actually, in [6] instead of varying the pressure at fixed chemical potentials \Xi 
and the temperature T, we varied the energy density under the condition that 
one should not vary it with respect to the particle occupation numbers. When 
E is varied, one should fix the particle densities rij and the entropy densities Si, 
which is equivalent to fixed particle occupations in our quasiparticle approach. 
Two procedures mentioned are equivalent, provided baryon-baryon hole and 
baryon-antibaryon excitation effects (loop contributions) are disregarded (as 
we did in (I12p - (j!5p ). An attempt to incorporate the baryon loop corrections 
into our scheme has been done in Appendix B. 

Now let us consider partial contributions to the pressure in eq. ffTTT) . 



3.2 The baryon contribution 



The contribution of the given baryon (antibaryon) species b G {b} to the 
pressure is as follows: 



f ^-P P^ ff) (7) 

P b [f, uq, R ] = -(2s 6 + 1) / -z— 2 V n b fi ch 



2ir 2 uj b 

u 

p=\0\, uJ b = \Jmf(f) +p 2 . (16) 

The spin factor s b = 1/2 for nucleons (N) and hyperons, while s& = 3/2 for 
A-resonances. 

The baryon set {b} to be used (taken from Table 1 of [Hj) was fixed above, see 
after eq. ([1]) . Little differences in masses of charged and neutral particles of the 
given species are ignored. Also we ignore small inhomogeneous Coulomb field 
effects and put the electric potential V = 0. The charge chemical potential 
/i c h is then related to the isospin composition of the system. For the isospin- 
symmetric system, N = Z, one has yU c h = 0. 

The Fermi-particle (baryon/antibaryon) occupation 



fb exp[K-K)/T] + l 
depends on the gauge-shifted values of the chemical potentials 
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A*6 — tb /^bar + £fc /4tr + ^ch ~ Qwb X^O ~ t b g p b Xp ^0 ■ 



(18) 



The baryon/antibaryon chemical potential of the fc-species is ^ = t& /ibar, and 
the corresponding strangeness term is /i^ = /i str . Baryon quantum numbers 
U, t b s , tf and are baryon charge, strangeness, isospin projection and electric 
charge, respectively, and proper charge conjugated values for antiparticles are 
given in Table 1 in [6]. 

3.3 Mean-field contribution 

It is convenient to introduce the coupling ratios 

Xmb = gmb/gmN, TH G {l7l} = 0~, LO, p, (19) 

and, instead of Xm, other variables 



since the pressure depends namely on this sort of combinations rather than 
on $ m and Xm separately. 

In terms of these new variables the contribution of mean fields to the pressure 
is as follows: 



Vm(f) = <U)/XIU) , 



(20) 



PTlf] 




Va(f) ~ U(f) 



(21) 



2CI 




(22) 



P p MF [/,^o] 




[dpN Xp Ro] ■ 



(23) 



Here the renormalized constants are 



a 



™>N g-mN 



(24) 



m 
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The net baryon density is given by [B|: 
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n B = , n b = (2s b + 1) / — — f b , (25) 

be{b} J Z7r 

where nf, is the baryon (antibaryon) number density and occupation baryon 
(antibaryon) density is defined by eq. (1171) . On the other hand, for fixed baryon 
species the contribution to the baryon density should obey the thermodynamic 
consistency condition 



dP 
dfi* b 



(26) 



Both quantities presented by fl26l) and fl25l) coincide provided contributions of 
boson excitations do not depend on the baryon loop terms (see Appendix B). 
In this case thermodynamic consistency of the model is preserved, see below 
in more detail. 



The isotopic charge density in the baryon sector is given by 



4=2E *6 n b x P b • (27) 
6e{6} 

The isovector baryon density n B plays the role of the source for the p-meson 
field = Ro . Therefore, for the iso-symmetrical matter (N = Z) one has 
n l B = and P p MF = . 

The net strangeness density of baryons and mesons reads 



n str = ]T tl m . (28) 

ie{6},{cx} 

Bearing in mind applications of the model to high-energy heavy ion collisions 
from AGS to RHIC energies we assume that all strange particles are trapped 
inside the fireball till the freeze-out. Therefore, the total strangeness is zero. 
Thus, we put n str = 0. This condition determines the value of the strangeness 
chemical potential fi str . 

Similarly, we may introduce the electric charge density 



^ch 



E 

i€{b},{ex} 



Q 



II; 



(29) 



The quantity n c h = {Z/A)ub determines the value of the charged chemical 
potential /i c h. For the symmetric matter, N = Z, ignoring Coulomb effects 
one may put fi n = fi p and /i c h = 0. 
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Our SHMC model pressure functional depends on four particular combina- 
tions of functions, r] a)P)UJ (f) and U(f). Note that the dependence on the scaling 
function r] a can always be presented as part of the new potential U obtained 
by means of the replacement U — > U + m 2 N c l (1 — 7)<r{f)) , and vice versa, so 
the potential U can be absorbed in the new quantity rj a . Thus actually only 
three independent functions enter into the pressure functional. Equation (fTTT) 
together with eqs. ffTSj) . f[2"Tj) . fl2"2"j) . f[2"3"j) demonstrates explicitly the equiva- 
lence of mean-field Lagrangians for constant fields with various parameters if 
they correspond to the same functions 7] P)U (f) and r] a (either U(f)) with the 
field / related to the scalar field a through eq. ([3]). In [6], we assumed r\ a = 1. 
Here we accept the same choice. 



3.4 Bosomc excitations 



To find the total pressure ffTTj) . one should define the contribution of bosonic 
excitations. Within our model and in agreement with [6] it is the sum of partial 
contributions 



Pbos.cx[/, ou , Rq, T] = P CT P art + pP art + P p P art + pP art 

ppart , ppart , ppart . ppart . ppart 
' r K + r r) "+* r K* "+" r i)> "+" ^ 

The pressure of the pion gas is 



(30) 



PP art = P.+ + P.O + P„- 



dp p 4 
2n 2 



x 



^7T+ (?) 



UJ n o(p) 



(31) 



Due to the absence of the u; — > 27r decay the coupling g*^ = 0. For N = Z, 
the field Po = and dependence of pion spectra on g* n disappears. Also, as 
in [6], we suppose g* n = ignoring a small pion mass shift. Then for both 
charged and neutral pions we may use 



u«±(p) = ^o(p) = yjml +p 2 . 
The pressure of the kaon gas is given as follows 



(32) 



ppart _ p 
1 

~ 3 



K+ + P K o + P K - + P R0 

00 



(33) 



dp p 
2tt 2 



f K +(uJ K +(p)) f K o(u K o(p)) 



uj k + (p) 
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fK-(VK-(p)) fK°(uK<>(p)) 

u K -{p) uro(p) 
where 



wjf±(p) = ±g* K uj ± g* K Ro + y]m*$ +p 2 , m* K = m K - g* K a (34) 

and g* mK = g m KXm for the parameters of the set B [6], we use here. Note that 
in neglecting Coulomb effects the energy of K + coincides with that for K° 
and the energy of K~ coincides with the K° energy. 

The ^-contribution to the pressure is given by 



+ 



dpp 4 



ppart _ 

n -3 



dpp 4 fr,(uT,(p)) 
2vr 2 u v (p) 



m 



*2 



(35) 



with 



,*2 2 



^ r lb /vf, iTr \ I / / 1 , K Vb 



- E ^ / 1+ E (36) 

6e{6} ■'tt / \ 6e{fe} •'t 



expressed in ref. [29] in terms of the total baryon scalar density J2be{b} n sb 



In the above equations the Bose distributions of excitations are 



/i= ! % 1 » ( 3? ) 

expK^m^ + p 2 -^)/^-! 

= + QiMch - gr ^ - Qrs^o, 

z G {bos.ex} = a',u', p' + , p'°, p'~; ir + , ir°, tt~; K + , K°, K~ , K°; 77; 

Here \ii = /i str for strange particles K and K* and /ij = — /i str for their 
anti-particles; Qi is the boson electric charge in proton charge units (ignoring 
Coulomb effects we put Qi = 0), QJ CC = +1 for particle, QJ ec = — 1 for an- 
tiparticle and QJ CC = for "neutral" particles (with all zero charges including 
strangeness). We take couplings as in [6], g^ = for all i except i = K,r). 
For 77 we select the parameter choice = 140 MeV, K v b = 0.2. As in [6] we 
assume m* K , = rriK*, m*, = rr?y, m* = m v due to the absence of the corre- 
sponding experimental data. In this paper, we will focus on the consideration 
of the isospin-symmetrical matter; therefore, we may put Rq = 0. 
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The density of the gas of Bose excitations of the given species i is determined 
by the integral 



9% J -y^-/*(p) > 1 e {bos.ex}, (3* 
o 



where gi is the corresponding degeneracy factor. 

Note that the p-wave pion and kaon terms can easily be included. For that one 
needs to replace u) n (p) and u>k(p) with more complicated expressions, see refs. 
[24"f26] . In order to obtain the temperature-density dependent pion and kaon 
spectral functions one needs to calculate the pion and the kaon self-energies in- 
cluding baryon-baryon hole loops and baryon-baryon correlation effects. These 
effects may result in appearance of the pion and antikaon condensates in dense 
and not too hot nuclear matter. Within the SHMC model we incorporate only 
interactions of particles through mean fields. Even in this case many coupling 
constants are not well fixed due to the lack of experimental data. Thus, we 
postpone the inclusion of the p-wave interactions to the future work. 

A more nontrivial task is to fix the terms pp art ; PP art , pp art . For the a' , 10 , p 
contributions to the pressure we use the standard ideal gas expressions with 
effective masses determined as follows. Let us first focus on the contribution 
of the a -excitations. In order to get pp art ; we should expand total pressure 
P[a, cuo(cr), Rq(ct)] in a = a — cr cl . The term linear in a does not give a 
contribution due to subsequent requirement of the pressure minimum in o~ cl . 
The quadratic term produces effective a particle mass squared, 



2 _ d 2 P[a,u; (a),R (a)] _ d"P[f,u (f),R (f)} ( df\ 
{m ° ] ~ do^ " dP \fo)' K } 

The first-order derivative dP/df = 0, as it follows from the equations of mo- 
tion. Since we deal with the strong interaction problem and the general solu- 
tion is impossible, one should use some approximations. First, keeping only 
quadratic terms in all thermodynamical quantities in boson fluctuating fields 
we disregard boson excitation contributions to the a , uj , p effective masses 
(higher order effects in boson fluctuating fields). Within this approximation 
we neglect the Pbos.ex. term in (1591) . Moreover, our aim is the construction 
of a thermodynamically consistent model where the baryon density and the 
entropy density, calculated by using the thermodynamic relation (126]) and re- 
lation 



dp 
b ~dr 



■ 

i' 



(40) 
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coincide with quasiparticle expressions (1251) and 



S = S hax + Sbos.cx., Sbar = J2 9b 7T 2 K 1 ~ .fiOM 1 ~ /&) - fb^fb], 

be{b} J Z7r 

Sbo,ex. = Y, 9ij [(1 + /<M1 + /*) - /an/,], (41) 

iGjbos.ex.} 

namely, the latter expressions in ref. [6]. In order to keep thermodynamic 
consistency of the model, we suppress the baryon contribution in This 
term is associated with taking into account of the baryon-antibaryon and 
baryon-baryon hole loops. Problems arisen, if one includes baryon loops, are 
discussed in Appendix B. 

Thus, the squared effective mass of the a' excitation will be given by the 
expression 



part* \ 2 



E 

m£{m) 



d 2 P m MF [/,. (/),fi (/)] (df 
df 2 [da 



(42) 



Using the relation x CT = $ 



/ we find that 



|_^v (1 _ /)2 (43) 
da m N y J J 

For the isospin-symmetrical nuclear matter N = Z, we obtain 



[m 



part* \ 2 



III 



N 



U ' + Cl 



(^) 2 (i-/r 

\m N J 



(44) 



One could use a different approximation introducing the effective mass term 
with the help of the expression for the Hamiltonian (Lagrangian). In the lat- 
ter case, baryon terms are added to the meson mean field contribution and 
derivatives of the Hamiltonian are taken at fixed ^ b- In this case 




K ffl t = (0) = ( ^ ) ( 5: ) . («) 

where H is the Hamiltonian and the averaging procedure is carried out over a 
thermal equilibrium state after taking the derivatives. We used that (^j^j = 0, 
as it follows from the equations of motion. Neglecting the meson excitation 
terms in H we find 



1(3 



6 



where the baryon (antibaryon) scalar density n s b is given by 



oo 9 

„, t = ( 2 St + i)/to/, 



Since in ref. [6] the value 7/ J 1 was chosen as a linear function of /, 



r]- 1 = (1 + z/)/(l + z/(n )), z = const, (47) 

(z = 0.65 in the KVOR-based SHMC model), we get £ be{b} ^ h ^U),t\ = 0j 
and both results 042p and (]45[) coincide. From here it is seen that in our 
approximation ( m P art *) 2 does not explicitly depend on the baryon variables 
and the temperature. Obviously, within this approximation our model remains 
thermodynamically consistent since meson excitation terms in the pressure do 
not contribute to the derivatives of the pressure over the baryon chemical 
potential (to obtain n B ) and over the temperature (to get entropy S). 

Note that in [5] we introduced the squared effective mass of the cx'-excitation 
as the second derivative of the energy density. Since at the same time the con- 
tributions of all baryon particle-particle hole and particle-antiparticle (loop) 
terms were suppressed, this result coincides with expressions and (Jj5j) . 

Within our approximation the masses of vector particles are given by 



[m 



( m P 



part*N2 / d H \ 9 2 P m [/, Up] , 

^ _ /d 2 H\ ^ g^L^g] _ n -I o q 



In this case zero-components of one u and three p-excitations prove to be the 
same as those following from the mean-field mass terms 



= m w |$ w (/)|, m£ art * = ro p |$ w (/)| . (49) 

Moreover, there are excitations of two magnetic-like components of uo and six 
magnetic-like components of p. Their masses are given by eqs. ( 1491) since there 
are no u and R mean fields, namely, expressions (14"9~]1 were used in [6]. 



17 



Thus, neglecting the baryon and boson excitation contributions to the effec- 
tive meson excitation masses we obtain a rather simple thermodynamically 
consistent description. 

3.5 Some numerical results 

In Figs. El [3] we demonstrate the size of corrections which arise provided boson 
excitations are incorporated in the equations of motion. Particle excitation 
masses are calculated following eq. (]42p for a meson excitation and (I49p for 
the nucleon, u>, and p excitations. Figure [2] demonstrates the ratio of the 
effective masses to the bare masses for the nucleon-u; -p (left panel) and for 
a (right panel) as a function of the temperature at n% = (thin curves) and 
at ub = 5 n (bold lines). Solid curves are calculations of the given paper, 
when boson excitations contribute to the equations of motion, whereas dashed 
curves show perturbative calculations of ref. [6j. Although for T > 120 MeV 
there appear pronounced differences between both treatments of excitations, 
the qualitative behavior remains unchanged. Effective masses of all excitations 
exhibit a similar behavior as functions of the temperature and the density, in 
a line with the mass-scaling hypothesis that we have exploited for the mean 
fields. 

Let us also compare our result for N, u, p effective masses (Fig. [2] left) with 
that previously obtained in the model of ref. [13], their Fig. 1 left (for N) and 
right (for uj, p). Shapes of the curves look similar. However in our case N, u, p 
are scaled by one scaling function, whereas ref. [13] used different scalings for 
N and for u,p. In their case the masses drop much stronger with the density 
increase but remain finite at high temperatures (T > 200 MeV), whereas in 
our case they drop to zero at T c ~ 210 MeV. 

The effect of the mentioned nonperturbative treatment of boson excitations on 
thermodynamic quantities is presented in Fig. [3j The temperature dependence 
is shown for the reduced pressure calculated in the given work (solid lines) and 
in the perturbative treatment of boson excitations [B] at ub = and ub = 5uq. 
One can see that the differences are rather noticeable only in the temperature 
interval 170 < T < 200 MeV, provided the baryon-meson couplings are not 
artificially suppressed. In the case (rig = 0), when all g a b couplings except for 
nucleons are artificially suppressed by factors 2/3 and 1/10, 0the temperature 
interval, where solid and dash curves deviate from each other, is broader but 
the value of the deviation is smaller. The curve for g ab suppressed by 1/10 
reasonably matches the lattice data for T c < T < 230 MeV. As has been 
mentioned in [6J one could fit the lattice data even in a broader region of 

1 As in Fig. 10 of [BJ, we suppress g^-couplings, not g m b, as mistakenly indicated 
there. 
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T, MeV T, MeV 

Fig. 2. Temperature dependence of the effective-to-bare mass ratio for nucleon-w'-p' 
(left panel) and for a (right panel) for ub = and ns = 5 tiq. The solid lines are 
our results, the dashed ones are perturbative calculations from [6]. To guide the eye, 
the horizontal dots show the m* = 2m n and m* = m n thresholds. 



Q_ 




T, MeV 



Fig. 3. Temperature dependence of the reduced pressure. Solid lines show cal- 
culations of the given work, whereas dashed lines demonstrate the results of the 
perturbative treatment of boson excitations [6]. The curves labeled by 2/3 and 1/10 
correspond to the case when all g a b couplings except for nucleons are suppressed by 
factors 2/3 and 1/10, respectively. Filled squares show the lattice QCD result for 
the 2+1 flavor case [30] . 

temperatures (up to 500 MeV) if one introduced Xa- < A violation of the 
universality of the a scaling would be also in a line with that we have used for 
uj and p, rju ^ 1 and r\ p ^ 1, as required to describe properties of neutron stars, 
cf. [26]. However we will not elaborate such a possibility in the present work. 
Instead we simply suppress g a b, b ^ N . Hydrodynamic simulations of heavy 
ion collisions are fully determined by the EoS to be described with quark and 
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gluon degrees of freedom or with only hadron ones. With such an EoS (with 
suppressed g^) a quark liquid would masquerade as a hadron one. 

180 
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^ 100 

80 
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40 

500 600 700 800 

M<g> MeV 

Fig. 4. Isentropic trajectories for central Au+Au collisions at different bombarding 
energies calculated in the present paper (solid lines) and with the perturbative 
treatment of boson excitations [6j (dashed lines). The results are presented for the 
SIS-to-AGS energies. Experimental points with error bars are taken from [31]. The 
freeze-out points marked by stars are obtained in [6j. Thin line corresponds to the 
freeze-out curve in [31] while dash-dotted line is that from [32]. 

In Fig. [4] we show the isentropic trajectories for central Au+Au collisions at 
different bombarding energies calculated in the present paper (solid lines) and 
with the perturbative treatment of boson excitations from ref. [H] (dashed 
lines). As we can see from Fig. [H for AGS energies the trajectories calcu- 
lated here almost coincide with those computed in ref. [6] where boson exci- 
tations were treated perturbatively. Then, as in [6] we extend our analysis to 
higher bombarding energies. The entropy per baryon participants was calcu- 
lated in [35] within the 3-fluid hydrodynamic model assuming occurrence of 
the first order phase transition to a quark-gluon plasma. The energy range 
from AGS to SPS was covered there. We use the S/N B values for E[ ab =158, 
80, 40 and 20 AGeV, at which the particle ratios were measured by the NA49 
collaboration (cf. [31]). At the RHIC we put S/Nb = 300 in accordance with 
an estimate in [3"4"] . 
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Fig. 5. Isentropic trajectories for central Au+Au collisions at different bombard- 
ing energies calculated in the present paper. Calculations are performed from AGS 
to RHIC energies for not suppressed couplings (solid lines) and with suppressed 
g a b couplings (except for nucleons) by a factor of 1/10 (dashed lines). The filled 
diamonds and triangles are obtained from the 47T particle ratios in [3TJ and [33], 
respectively. Filled circles are RHIC data based on the middle rapidity particle ra- 
tios |31] . Open circles, squares and triangles are the lattice 2-flavor QCD results [M] 
for S/Nb = 30, 45 and 300, respectively. The thin line corresponds to the freeze-out 
curve in [3T] . 

Since there is a reasonable fit of the lattice data on the pressure for T > 
200 MeV with reduced meson-baryon coupling constants, see [6] and Fig. [3j 
for high collision energies we show in Fig. [5] both the results, when couplings 
are not suppressed (solid curves) and also when all g a b couplings except for 
nucleons are suppressed by a factor of 1/10. As is seen, the decrease in 
improves the high-temperature description (as compared to the case of not 
suppressed couplings) for the lattice data at high temperatures not only for 
the case //& = but also for fi b ^ 0. For low temperatures the solid and 
dashed curves (the curves for not suppressed and suppressed g a b couplings, 
respectively) are very close to each other. Also note that one should be rather 
critical performing comparison with the existing lattice data, especially for 
the case fif, ^ (presented in Fig. [5] for 2 flavors). Doing this comparison we 
tentatively hope that results of future more realistic lattice calculations will 
not deviate much from the existing ones. 
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In all cases our new calculations presented in Fig. [5] (if one additionally sup- 
presses couplings), as well as corresponding perturbative calculations [6], de- 
scribe reasonably the lattice data on the T — ps trajectories and freeze out 
points. 

Concluding, our improvement of the model with keeping boson excitation 
terms in the equations of motion neither spoils nor improves the agreement 
with the lattice results and with the thermodynamic parameters extracted at 
freeze out. Thus, both perturbative and nonperturbative treatments of bo- 
son excitation effects can be used with equal success. However, one should 
note that thermodynamic consistency conditions are fulfilled exactly in the 
non-perturbative treatment of the given work, whereas in the perturbative 
treatment of [B] they were satisfied only approximately. 



4 Inclusion of finite resonance widths 

In the framework of our model we treated all resonances as quasiparticles 
neglecting their widths. Excitations interact only with the mean fields. Thus 
imaginary parts of the self-energies as well as some contributions to the real 
parts are not taken into account. Definitely it is only a rough approximation 
carried out just for simplification. 

At the resonance peak the vacuum A-isobar mass width is r^ ax ~ 115 MeV. 
In reality Ia is the temperature-, density- and energy-momentum-dependent 
quantity. For low A-energies the width is much less than T^ ax . A typical A 
energy is loa — m* A ~ T. Thus for low temperature, T < ep (e? is the nucleon 
Fermi energy) the effective value of the A-width is significantly less than 
-pmax these temperatures the quasiparticle approximation does not work 
for A's but their contribution to thermodynamic quantities is small. When 
the temperature is > m n there appears essential temperature contribution to 
the width and the resonance becomes broader [36J. At such temperatures A's 
essentially contribute to thermodynamic quantities. Only for temperatures T > 
T^ ax (T) the quasiparticle approximation becomes a reasonable approximation. 

The p- and cr-mesons also have rather broad widths. E.g., at a maximum the 
p-meson width is about 150 MeV. The observed enhancement of the dilepton 
production at CERN, in particular in the recent NA60 experiment [37] on 
p + p~ production, can be explained by significant broadening of the p in matter 
[2], though decreasing of the p mass could also help in explanation of the data 
[38_[H Besides, as was shown in [30] to be consistent with the QCD sum rules, 

2 As demonstrated in [38], the calculated large mass shift is mainly caused by the 
assumed temperature dependence of the in-medium mass. Inclusion of this tem- 
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both the collisional broadening and dropping of the p mass should be taken 
into account. Even more, one should be careful with interpretation of the 
NA60 experiment: Dileptons carry direct information on the p meson spectral 
function only if the vector dominance is valid but generally it is not the case 

m- 

Also particles which have no widths in vacuum like nucleons acquire the widths 
in matter due to collisional broadening. Their widths grow with the temper- 
ature increase, cf. [12]. As we have argued above in case with A isobars, the 
quasiparticle approximation may become a reasonable approximation at suf- 
ficiently high temperature, if T > T(T). 

Let us roughly estimate the effect of finite baryon and meson resonance widths 
on particle distributions and the EoS in order to understand how much and in 
which temperature-density regions these effects may affect results calculated 
within our quasiparticle SHMC model. 

It is convenient to introduce single-particle spectral A r and width T r functions 
(operators), cf. [12], 



A x = -2ImGf (q) = -21m ^ 1 ^ , f r = -2ImSf , (50) 

r v ; M r + ir r /2 r v ' 

where G R (q) is the full retarded Green function of the fermion and Y> R is the 
retarded self-energy. The quantity 



M r = (G^ R Y l - ReSf (51) 

demonstrates the deviation from the mass shell: M r = on the quasiparticle 
mass shell in the matter, G 0,R is the free Green function. Following defini- 
tion (1501) the spectral function of the fermion has dimensionality of m~ l and 
the width, dimensionality m, whereas for bosons the spectral function has 
dimensionality of m~ 2 and the width, dimensionality m 2 . 

Let us start with the consideration of a fermion spin 1/2 resonance (superscript 
f). The spectral function satisfies the sum rule: 



-Tr 



7o 



Al (+) (q Q ,q) 



^r(-)0o, 



dqo 
2tt 



(52) 



perature dependence modifies the scaling hypothesis originally claimed by Brown 
and Rho. Some arguments on what the proper mass-scaling predicts for dilepton 
production in HIC, e.g. NA60, were given in |39j. 
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7o is the corresponding Dirac matrix, and subscripts (±) specify particle and 
antiparticle terms. The trace is taken over spin degrees of freedom. 

Simplifying the spin structure, as it is seen from (150]) . we can present 



-Tr 7o ^ (+) = 4 = (53) 

Dealing further with A f T — A T iq we may not care anymore about the spin 
structure. The A spin 3/2 resonance can be considered in a similar way, see 

m- 

It was argued in refs. [B] that at a self-consistent description the particle 
densities are given by the Noether quantities: 



v ' 

Here fi* = fi* N for a nucleon resonance such as A. The antiparticle density is 
obtained from (154|) with the help of the replacement fi* — > —fi*. The baryon 
density of the given species is then n]? ar = n T (fi*) — n T (—fi*). N r is the degen- 
eracy factor. 

To do the problem tractable, instead of solving a complete set of the Dyson 
equations, we may select a simplified phenomenological expression for A\ (com- 
pare with [13]), e.g., 



Af 2£u \2YJs) + 26} , 9 

Al = 1 rV L — l - , f = const, s = u 2 -p 2 >0, 55 

r (s - mff + \T T {s) + 5} 2 

with 5 — > +0. The value 5 — > is introduced to easier extract the quasiparticle 
term. Separating out the quasiparticle pole we obtain 



4 = , .2 2 ; l ~ J , , 12 + 2 ^ • 2n5(s - K) 2 ) «(«u, " K) 2 ). (56) 

Here s t h is the resonance threshold value of s. Note that within our ansatz, 
the spectral function depends only on the s-variable. It might be the case 
only for a dilute matter, when the density and temperature dependence of 
the width is rather weak. Furthermore, instead of a calculation of the density- 
temperature dependent part of the width, which actually can't be performed in 
the framework of our model, we vary the energy dependence and the amplitude 
of the width thus simulating in such a way collision broadening effects. 
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For the decay of the resonance into two particles (r — > 1 + 2) one may use a 
simple s- variable dependence of the width: 



r, W = r. m, F (s) ( ^Xt)) "^ - K + m ^ • (57) 

n 2 ( - * ~ _ (s - (ml + m* 2 ) 2 ){s - (ml - m*) 2 ) 

Here T = const, the width tends to zero at the threshold s — > s t h = (tj* + 
mj) 2 , a = / + 1/2, with a = 1/2 for s and a =3/2 for p resonance. An extra 
form-factor, F(s), is introduced to correct the high-energy behavior of the 
width. 

To simplify expression fl57|) . we may expand r r (s) near the threshold trans- 
porting remaining s-dependence to the form factor: 



r r (s) = T F(s) m T 



f 5 th 

1/2 



0(s - s t h) 



(58) 



where one can take 



F 1 + [(s - s th )/s )P (59) 

with so and (3 being constants. The parameters can be fitted to satisfy exper- 
imental data. 



The energy dependence of the width causes a problem. With a simple ansatz 
for the behavior T r (s) we get a complicated m* 2 (s) dependence, as it follows 
from the Kramers- Kronig relation. However, since m* 2 (s) is a smooth function 
of s, one may ignore this complexity taking for simplicity 777/ j. £lS Si constant. 
The factor £ is introduced to fulfill the sum-rule: 



/J£ = l, (60) 



that yields £ ~ 1 + O(r /m*) (for m* > T ), m* > m* + m^. In the case 
of m* < m* + Trig there appears an extra quasiparticle term in the spectral 
function, that contributes to the sum-rule. 

In Figs. [6]and[7J we present the ratio Ra = n r £ s /n A p of the A-isobar density, 
calculated following eq. (!54{) . to the quasiparticle density, as a function of the 
temperature. The results are presented for the baryon density % = (in Fig. 
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Fig. 6. Ratio of the A isobar density, calculated with the inclusion of the width, 
to that of the quasiparticle one, = rf£ s /n A p , as a function of temperature at 

1/2 

us = for two values of the resonance threshold energy s t ^ . The parameters of 
calculation are presented in the figure. 



[6]) and for % = 5no (in Fig. [7]). Our aim here is to demonstrate the effect of 
a finite particle width. Therefore, instead of searching for the best fit of the 
spectral function to available experimental data we vary the parameters to 
show how strongly the density ratio may depend on them. We take into account 
the p-wave nature of the resonance and use Tq = 115 MeV. The threshold 
quantities are chosen to be s\[* = m* N (left panels) and s][ 2 = m* N + m n (right 
panels). In vacuum = m N + m n . Thus, taking = m^+m^ we simulate 
the vacuum resonance placed in the mean field (in our model m* = m^). With 
s\(^ = m* N we simulate the effect of in-medium off-shell pions (virtual pions 
can be produced in matter at any energy). The form factor F is computed 
with /3 = 3 and s = (500 MeV) 2 (solid lines) and s = (1000 MeV) 2 (dash 
lines) to present the dependence of i?A on the high-energy behavior of the 
width, which is not well defined even in vacuum. 

As we can see, in all examples the curves are rather flat in the temperature 
range T > 50 MeV-fdOO MeV. For n B = 0, at T ~ 170 MeV there appears 
a slight bend associated with a sharp decrease in the nucleon effective mass 
for T > 170 MeV. For % = 5 n the bend is smeared. A substantial devia- 
tion of the Ra ratio from unity for T > 50 MeV-4-100 MeV in the case with 
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Fig. 7. The same as in Fig. [6] but for n# = 5 uq. 



the cut-off So = (1000 MeV) 2 is due to a larger width at high resonance en- 
ergies (for large s) compared to the case with the cut-off so = (500 MeV) 2 . 
In the latter example, the quasiparticle approximation becomes appropriate 
for T > 100 MeV. On the contrary, for s = (1000 MeV) 2 and n B = the 
quasiparticle approximation does not work at all. For low temperatures, the 
i?A ratio is significantly higher than unity and for higher temperatures it be- 
comes essentially smaller than unity. The broader the width distribution is, 
the smaller i?A at large temperatures. The ratio dependence on the thresh- 
old value s]^ 2 = m* N is rather pronounced for s = (1000 MeV) 2 but it is 
only minor for s — (500 MeV) 2 . Thus, we conclude that taking into account 
the energy dependence of the width might be quite important for very broad 
resonances, when the width only slowly decreases with energy. If the width 
drops rather rapidly with the energy increase, the quasiparticle approximation 
becomes appropriate for calculation of thermodynamic quantities already at 
not too high temperature. The baryon density dependence of the ratio i?A is 
not so pronounced (especially for s = (500 MeV) 2 ), since in our parameter- 
ization the spectral function depends on the density only through the value 
m*(riB) and the choice of s t h- For low temperatures (T < 50 -j- 100 MeV), the 
i?A ratio becomes significantly larger than unity. We also pay attention to 
the shift of the reaction thresholds due to the dependence of the width on s 
and the threshold s t h on the density and temperature. This point can be very 
important for fitting of particle momentum distributions. 
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To demonstrate the effect of the finite resonance width on thermodynamic 
characteristics of the system, we calculate the energy density of the non- 
interacting resonances (however with width). Then 



oo 



A-irp 2 dp r ds 



K = N r 



J 



(2tt) 3 J 4vr 



uA T f T + (/i r 



/*r), 



(61) 







where the degeneracy factor for A is iV r = 16. 

In Figs. [8] and [9], we show the ratios of the energy density for A's (with and 
without width) to the total baryon energy density at ns = and ns = 5 no, 
respectively. The solid and dashed curves correspond to calculational results 
for the width with s = (500 MeV) 2 and s = (1 GeV) 2 . The dash-dotted 
curve is computed within the quasiparticle approximation. Calculations are 

1 /2 

performed for two values of the threshold energies s t ^ = m* N (left panels) and 
sl[ 2 = m* N + m n (right panels) . 

As was expected, in all cases the quasiparticle result is much closer to that 
for sq = (500 MeV) 2 than for sq — (1 GeV) 2 . In the former case, the differ- 
ences are almost negligible. For s = (1 GeV) 2 the ratio remains smaller than 
the quasiparticle result, except for low temperatures. Differences in the ratios 
with and without taking into account the width in all cases are not too no- 
ticeable. The density dependence of the ratio proves to be pronounced even in 
our model, although the density dependence of the width is not incorporated 
explicitly. Summarizing, for calculation of thermodynamic characteristics, one 
may use the quasiparticle approximation for baryons in the whole temperature 
interval under consideration provided the high energetic tail of the resonance 
width-function is not too long. If a resonance has a long energetic tail, then 
the longer the width tail is, the more suppressed the ratio of the resonance 
energy to the total energy, as compared to the corresponding quasiparticle 
ratio. 

For charged boson£i| the spectral function follows the sum- rule, cf. |45j . 



We again consider a dilute matter assuming that the spectral function depends 
only on the s-variable. Only in this case one may consider a single spectral 
function for vector mesons, like u and p, whereas in the general case one should 
introduce transversal and longitudinal components. 

3 As before, by the charge we mean any conserved quantity like electric charge, 
strangeness, etc. 




(62) 
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Fig. 8. Ratio of the A isobar energy density, calculated with the inclusion of the 
width (-Ea)j to the total baryon energy density (Eg), as a function of temperature 

1/2 

&t ub = for two values of the resonance threshold energy s./ . The parameters 
of calculation are presented in the figure. E q £ is calculated within the quasiparticle 
approximation . 

The Noether charged boson density (of given charge) is given by 



o v ; 

For practical calculations it is convenient to present the spectral function in 
the form: 



^ = t7^4mm + f ' 2rf(s " K)2) " (s,h " K)2) (64) 

with £ = const introduced to fulfill the sum- rule. Replacing T r (s) = |r]?(s), 
we may use eq. (157|) for T T (s) with a = 1/2 for s and a = 3/2 for the p 
resonance. 

In Figs. [10] and [Til the ratio R p = n p /n^ p of the p + density, calculated following 
eq. ([531 . to the quasiparticle density is shown as a function of the temperature. 
Here we take into account the p-wave nature of the resonance and use two 
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T, MeV T, MeV 

Fig. 9. The same as in Fig. [H]but for ub = 5no- 
values for the width: a small width Tq = 75 MeV (left panels) and a very 

1 /2 

large one T = 250 MeV (right panels). We use s t £ = 2m 7r , thus taking 
ml = m* 2 = m n . The results are presented for ub = and ub = 5 no in Figs. 
[10] and [UJ, respectively. These figures show only a moderate dependence of the 
ratios R p on the value of the width at the resonance peak (on r ). At high 
energies the energy dependence of the width is more pronounced (compare 
solid and dash curves). In our model the density dependence of R p reflects 
the behavior of m*(riB)- For ns = the ratio R p becomes larger than unity 
for T < 100 MeV, whereas for ub = 5 no R p is larger than unity only for 
T < 20 4- 50 MeV. The ratio R p {ub = 5 n ) < R p (ub = 0) since the resonance 
mass m* p {n B ) for n B = 5 n is closer to the threshold value 2m n than in the 
case of Ub — 0. 

Slight bends of the curves at T ~ 170 MeV are associated with dropping of 
the effective mass of p below the threshold. Then there arises a quasiparticle 
contribution to the spectral function being added to the high-energy width 
term (second term in (1641) ). In general, the behavior of R p and i?A is simi- 
lar. The quasiparticle approximation is rather appropriate for T > 100 MeV, 
provided the widths have not too long high-energy tails (see curves with 

1 /2 

s = 500 MeV). Under these conditions our quasiparticle SHMC model works 
well. If the resonance width has a very broad high-energy tail (see curves with 

1 /2 

s = 1000 MeV), the deviation of the i? r ratio (for the r-resonance) from 
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Fig. 10. Ratio of the p meson density calculated with the inclusion of the width to 
the quasiparticle one, R p = n^ es /n^ p , as a function of the temperature at ns = 0. 
The parameters of calculation are presented in the figure. 

unity is rather pronounced, even for T > 100 MeV. Since in this case R T < 1, 
when broad resonances are included, there may appear a possibility to match 
the lattice data for T > 170 MeV even without suppressing the coupling con- 
stants. The latter procedure was used in [B] to demonstrate a possibility to fit 
the lattice results within our quasiparticle SHMC model. 

Broad hadronic resonances having a short lifetime are of a particular inter- 
est for dynamics at the late stage of relativistic heavy-ion collisions. As the 
system expands and cools, it will hadronize and chemically freeze out (van- 
ishing inelastic collisions, no creation of new particles). After some period of 
hadronic elastic interactions, the system reaches the kinetic freeze out stage, 
when all hadrons stop interacting at all (vanishing even elastic collisions). Af- 
ter the stage of the kinetic freeze out, particles overcoming remaining mean 
fields free-stream towards the detectors, where measurements are performed. 
The lifetimes of the p meson and A-baryon are r p ~ 1/Top — 1-3 fm/c 
and ta — 1.7 fm/c, respectively, being small with respect to the lifetime of 
the expanding system. So short-lived resonances can decay and regenerate in 
scattering process all the way through the kinetic freeze out. The regeneration 
process depends on the hadronic cross sections of resonance daughters. Thus 
the study of different resonances can provide an important probe of the time 
evolution of the source from the beginning of the chemical freeze out to the 
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Fig. 11. 


The same as in Fig 


[JO] but for 


n B = 5 n . 



kinetic one. In-medium effects can modify properties of hadron quasiparticles 
and resonances during this stage. In particular, it concerns finite widths and 
energy-momentum relation evolving in time. Values of chemical potentials and 
the temperature characterizing particle momentum distributions become com- 
pletely frozen at the kinetic freeze out. From this moment only the mean fields 
can survive and the hadron effective masses further evolve towards the bare 
masses with which particles reach detectors. 

In Fig.[12]the ratios of yields of resonances to the yields of stable particles with 
similar quark content, the p / 71 " - and A ++ /p ratios, are presented as a func- 
tion of temperature. For all four species considered, our calculations take into 
account the feed-down from higher resonance^ n{ eed = rii+Y^r n r T r ^ l /T r . Be- 
ing in an agreement with experiment, the statistical analysis of stable hadron 
ratios at the chemical freeze-out [21] shows that these ratios are getting almost 
energy-independent at ^snn - 100 GeV. This statement seems to be valid 
also for resonances provided they are treated within the quasiparticle approx- 
imation, as follows from the weak T dependence of the ideal gas (IG) model 



4 Those hadrons from the particle data table which originally were not included into 
the SHMC model set are treated here as the ideal gas of resonances with vacuum 
masses and vanishing widths. 

5 It is a quasiparticle model with the vacuum masses for all hadrons. Mesons with 
masses mi < 1.6 GeV and baryons with mi < 2.2 GeV are included. 



32 



0,6 



s Q =(500 MeV) 2 



CO 

o 

-t—> 
CO 



0,1 -j. 



0,5 - 



0,4 - 



0,3 - 



0,2 -• ► 
Tj. 



T 




0,0 



120 



130 



140 



150 



160 



170 



180 



T, MeV 



Fig. 12. The resonance ratios as a function of temperature. The solid lines are cal- 
culated within the SHMC model with accounting for the resonance width according 
to eqs. (|54p . (|58p . Parameters are the same as in Figs. [61 and 1101 (right panels, solid 
lines). The dashed lines are calculated for the IG model. These results are obtained 
at fig = 20 MeV for the RHIC energy. Experimental data [38] for central d + Au 
collisions are plotted by straight dotted lines with error bars. 

results (see dashed lines in Fig. [T2l) . As is seen in Fig. [T2], the resonance ratios 
calculated in the SHMC and IG models almost coincide with each other till 
temperature about 140 and 155 MeV for the A ++ jp and p°/vr~ ratios, respec- 
tively, and then the difference between them increases with T. The growth of 
the resonance yield in the SHMC model is mainly due to the dropping of the 
effective mass at high temperatures. Dependence of the resonance abundance 
on the value of the width and other parameters (at the same temperature) is 
rather moderate (compare the solid and dashed curves in Figs. ISl and fTOi) . 

Earlier, the resonance statistical treatment has been considered in refs. [46] 
and |47|. Both models are ideal gas ones but in the hadronization statistical 
model [17] the vacuum (energy independent) resonance widths are included. 
At the same values T, /ig these models provide quite close results but if one 
refers to the particular collision energy, their results differ due to making 
use of different approximations for the freeze out T and Hb as functions of 
the collision energy (see Fig. [5]). In the hadronization statistical model |47j . 
predictions were made only for the SPS energies. 

Very recently the short-lived resonances A ++ and p° have been measured by 
the STAR Collaboration at RHIC in d+Au collisions [32]. The observed p / 71 " 
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and A ++ jp ratios practically do not depend on the charged particle number at 
the mid rapidity, dN c h/dri (i.e. on centrality), and for central (20% centrality) 
collisions are 0.139 ± 0.050 and 0.206 ± 0.034, respectively [IB]. As is seen in 
Fig. O calculated results for the SHMC model cross the experimental line at 
T ~ 145 and ~ 160 MeV0for A ++ jp and p°/7r~, respectively. However, large 
experimental error bars do not allow to make preference to any specific model, 
though the IG model predictions are slightly but regularly below that for the 
SHMC model. One should remind that at the RHIC energy the chemical freeze 
out temperature is estimated as 177 MeV [4TIT45] and kinetic one is about 120 
MeV. The p°/vr~ ratio was also measured in peripheral Au + Au(200 GeV) 
collisions @H] to be as large as p°/n- = 0.169±0.003(stat) ±0.027(syst). Note 
that the existence of the difference between the experimental value of p°/tt~ 
ratio and the statistical model result has been considered as a problem in 
ref. [UJ. Following their estimate, p°/vr~ ratio is about 4 • 10 -4 at the kinetic 
freeze out temperature 120 MeV, and it is 0.11 at the chemical freeze out 
T = 177 MeV. Our result (~ 0.06 at T = 120 MeV) strongly deviates from 
the statistical estimate [17] at the kinetic freeze out. At the chemical freeze 
out temperature, T = 177 MeV, predictions of the p°/tt~ ratio practically 
coincide in all ideal gas based models remaining lower than the experimental 
ratio, whereas our SHMC model is able to reproduce experiment at T ~ 
165 -=- 170 MeV. It is also of interest that the measured resonances have a 
mass shift about 50-70 MeV [H] which is consistent with the SHMC model 
results at T ~ 160 MeV as presented in Fig. [2j Thus, one may infer that due 
to a possible decay, regeneration and rescattering, the broad resonances, if 
they are treated within the statistical picture with taken into account their 
in-medium mass shift, freeze at a somewhat lower temperature than that at 
the chemical freeze out but this temperature is significantly higher than the 
value of T ~ 120 MeV, characterizing the kinetic freeze out [50] . A dynamical 
consideration of the freeze out is needed to draw more definite conclusions. 



5 Concluding remarks 

In this paper we made attempts to find several improvements of the SHMC 
model of [6]. 

In [6] the boson excitation effects were assumed to be small. Therefore their 
contribution to thermodynamic characteristics was calculated using pertur- 
bation theory in the fields of boson excitations. Thereby, boson excitation 
contributions in the equations of motion were dropped. The approximation 

6 The temperature concept should be used with care for d + Au collisions. In exper- 
iments, temperature is usually associated with the inverse slope of transverse mass 
distributions. 
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made needs justification. Here we incorporated the boson excitation terms in 
the equations of motion and then calculated thermodynamic quantities (in- 
cluding boson excitation parts). Corrections to the effective masses of the a, 
and a;-p-nucleon excitations turn to be minor for T < 100 -j- 120 MeV and grow 
with the temperature increase. Corrections to thermodynamic characteristics, 
like total pressure, energy, entropy etc., remain moderate even at higher tem- 
peratures. Qualitatively, one may conclude that all results of [6] remained 
unchanged. With boson excitation terms incorporated into the equations of 
motion, our quasiparticle SHMC model fulfills exactly the thermodynamic 
consistency conditions. 

Then we discussed possible effects of the resonance widths. We assumed vac- 
uum but energy- dependent widths of resonances. Under this assumption the 
width effects included do not change qualitative behavior of the system. Nev- 
ertheless, one should note that in dense and/or hot matter particle widths 
may acquire essential density- and temperature-dependent contributions that 
may significantly affect properties of the system, e.g., see [12], where for the 
case of hot baryon-less system it is shown that width effects may completely 
smear fermion distributions. We illustrated that the estimated yields of short- 
lived resonances to be important for the late stage of relativistic heavy-ion 
collisions can be described sufficiently well if one treats those resonances in 
the framework of our SHMC model with the vacuum widths at the freeze out. 

In the paper, the a variable was considered as an order parameter. Then the 
effective masses of a , u , p excitations and the effective masses of nucleons 
have a similar behavior as a function of the baryon density and tempera- 
ture. The effective masses remain rather flat functions of the temperature and 
sharply drop to zero only in the vicinity of T c . In the large temperature inter- 
val mentioned, the effective masses as a function of the density first decrease 
with the density increase and then, at a very high density, begin to grow, cf. 

i- 

In Appendix A we calculated what modifications of our model could be, if 
the a-u and p fields were treated on equal footing. In the latter case, the 
density-temperature behavior of the a excitation mass essentially differs from 
that of the u -p excitation masses, namely, the a excitation mass drops to 
zero at n c ~ 4.2 and only slightly depends on T. Thus, if this model was 
applied to analyze heavy ion collisions, we would face with a problem of Bose 
condensation of a excitations for ub > n c ~ 4 4- 5 in a broad temperature 
interval. Another unpleasant feature of a model version like this is a drastic 
difference in the density behavior of the a effective mass and that of uj -p 
excitations. These are additional arguments in favor of the treatment of the 
a variable as an order parameter (as we did earlier in ref. [6] and here in the 
paper body). 



35 



Then in Appendix B an attempt was made to incorporate the nucleon-nucleon 
hole and nucleon-antinucleon loop effects in our model. The loop terms men- 
tioned were calculated within the perturbative approach. If these terms are 
taken into account, the a excitation effective mass exhibits rather unrealistic 
behavior. Thus, the model completely loses its attractiveness if baryon loops 
are included (at least within the perturbative approach). We have checked that 
this unpleasant feature is a common feature of many RMF models, including 
the original Walecka model (previously authors of ref. [51] arrived at a similar 
conclusion considering vacuum loop corrections in the Walecka model). In the 
Fermi liquid approach to diminish contribution of the fermion loop terms one 
incorporates the vertices corrected by a short-range baryon-baryon interaction 
introduced with the help of the Landau-Migdal parameters, cf. [24 25 26]. In 
the framework of our SHMC model short-range correlation effects are simu- 
lated by the ou, a, p exchanges and are non-local. Thereby, we do not include 
these effects in the present work. Summarizing, either higher order fluctuation 
effects should be included in all orders, that is a complicated problem, or they 
should be skipped within RMF based models. So we skipped these terms in 
our truncated scheme. 

Concluding, the SHMC model introduced in ref. [6] can be considered as a 
reasonable model for application to the description of hadronic matter in a 
broad baryon density-temperature range, provided higher order fluctuations 
of fermion fields are not included. 
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Appendix A. Treatment of the a field as an independent variable. As 

noted in the paper, we continue to suppress contributions of baryon loops. If 
the a and ljq fields are treated on equal footing, i.e. as independent variables 
("ind.var."), one can determine the effective a -mass using partial derivatives 
of the pressure, 



/ m part*\2 _ d 2 E m g{m} F [f , ^o] 

\ m a J ind.var. ~ ~Qp | 

rather than full derivatives. 
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It is convenient to present 



m.e{m} Pfn m.e{m} 



dP df 



2 



a- 



pjl pMF / rs , \ 2 /12 y> pMF n , 



^ 2 [df) +2 du df df {bb} 

Taking derivatives in ( 1651) and ( 1661) we use the exact equation of motion for / 
and uj , and suppress the boson excitation and baryon-loop contributions to 
the pressure in the final expression. 

Taking partial derivatives we get 

- P^wjTVM] ^ < + ^ _ ulml (67) 

and 5 a = — co^m 2 , where rj a = 1 was used. 

If the effective a'-excitation mass squared is determined using the Hamiltonian 
we obtain 

\ Qp J = ~cT + U f ~ u o m " + UG ~dp ^ 9ubtbnb, (68) 

that differs from (ISTl) by the last term. 

In FigfT3l we show the ratio of the effective-to-bare masses of the a' excitation 
as a function of the baryon density at different values of the temperature for 
two possible treatments of the a- field, as an order parameter, see (]42i) . and 
as an independent variable. In the latter case two expressions, ( 1651) and (168|) . 
were used to calculate the value m^ 1 *. The solid lines are evaluated following 
eq. (f65]) for T = 50, 150 and 175 MeV (from the top to the bottom); and the 
dashed curves, using (168j) . In both the cases the effective a excitation mass, 
calculated following (1651) . drops to zero at % » 4.2no (for T = 50 MeV) 
and for ub ~ 3.8 no (for T = 150 MeV), respectively, which demonstrates a 
moderate temperature dependence of the effective mass up to T ~ 150 MeV. 
For higher temperatures the critical density significantly decreases. In all the 
cases, the differences between calculations following 065)) and 068|) are mi- 
nor; however, they significantly deviate from calculational results following 
eq. (f42j) (dash-dotted curves). The behavior of the effective a excitation mass 
calculated using (165!) and (168!) is in contrast with the behavior of the to -p - 
N effective masses which do not reach zero for all densities (see dash-dotted 
curves). Since one of our main goals was to construct a model based on the 
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Fig. 13. The effective a excitation mass as a function of density calculated using 
eqs. (|65p (solid curves), (|68p (dash lines), and (|42p (dash-dotted curves) for different 
temperatures. Straight dotted lines correspond to m a = 27717,- and m a = m^. 

idea of a rather similar behavior for all o'-uJ-p-N excitation masses, we refuse 
considering the a field, as an independent variable, i.e., we refuse to consider 
it on equal footing with the uj and p fields. Therefore, instead of either eq. fl65l) 
or eq. fl68|) . here we use eq. (I42p treating the a field as an order parameter. 

Appendix B. Inclusion of baryon loops. 



In the general case, the a excitation mass is given by eq. (139]) . provided the a 
field is treated as an order parameter. We will continue to keep only quadratic 
terms in fluctuating boson fields. Thus, we drop the boson excitation term 
in the pressure in expression (13"§1) but we will keep the baryon term. So in 
contrast with (|4"2|) . we use the expression 



part* \ 2 



F [/^o(/)] 

dP 
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+ 



df 2 



'df 
da 



(69) 



An important difference between (|6T?|) and (141)1) (the latter expression yields 
the same result as (142]) ) is that derivatives of the Hamiltonian are taken at 
fixed ^b, whereas the total pressure P depends on the baryon occupations, 
which should be varied. Thereby, (1691) includes extra contributions from the 
baryon-baryon hole and baryon-antibaryon excitations. 

Let us now calculate an additional, purely baryon contribution to the a - 
excitation mass incorporated in (|69|) . Using that 



df du b 



df 



df 



df b df b 



p dp dujb 



and (|26|) . with the help of partial differentiations we find 



(70) 
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We also use the equations of motion and the relation (^r) 
Here the quantities 



(71) 



-x ab m N . 



B h 



L b H * ) 



and 



(72) 



L b = N b j i^^b fb (73) 
o 

are the baryon-baryon hole and baryon-antibaryon loop-terms taken at zero 
incoming energy and momentum. 

We also introduce a similar quantity 
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Fig. 14. The ratio of the effective a' excitation mass to the bare mass as a function 
of density (top) and temperature (bottom) calculated by means of eqs. (|4"2j) . (solid 
curves), eq. (|78j) (dashed curves) and eq. (fTTJ) (dash-dotted lines), see (|69|) . The 
straight dotted lines correspond to the level of m a = Im-^ and ra a = m^. 
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(74) 



Then we calculate an additional contribution 



d 2 P[f,uj (f)} , d 2 P[f,u ] 



dp 



dp 
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(75) 



that distinguishes full and partial derivative terms. Here 
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(76) 



with (rrC rt *-°f = d 2 P MF /dul 



In [6], the baryon excitation loop L b contributions were suppressed. On the 
other hand, one may check that L b > n sb /ml, L b > L b . Therefore, in [6] we 
disregarded the term 5^(mP art *) 2 as the whole, as well as the terms oc B b and 
L b in (1751) and (JT_6j) . Thus in [6] eq. (T4"2]) is actually used, being compatible 
with (T451). 

Incorporating the loop terms mentioned, one should use the following expres- 
sion 



part* \ 2 
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(77) 



where partial contributions are given by ( 1671) . ( 1711 and ( 1751) . 



If ex and u;o are treated on equal footing, i.e. as independent variables ("ind.var." ), 
one should use 



part* \ 2 



d 2 P[a,u; ] d 2 P[f,u ] (df 



'ind.var. 



da 2 



dp 



da 



(7S 



The ratio of effective-to-bare masses of the cx'-excitation is shown in Fig. [14] as 
a function of the baryon density for three values of temperatures (top panel) 
and as a function of the temperature for three values of the baryon density % 
(bottom panel) calculated by means of eq. (J77J), i.e., following (1691) . Results 
are compared with (178]) and (142]) . As we can see, the inclusion of the baryon 
loop terms (performed within the perturbation theory) completely destroys 
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all achievements of our SHMC model (besides the case % = when these 
loop terms are suppressed). 

Thus, at the moment we can see no simple way to generalize the mean field 
based the SHMC model which would include baryon-antibaryon and baryon 
- baryon hole fluctuations. 
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